Exact results for quantum phase transitions in random XY spin chains 
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The effect of disorder on the quantum phase transitions induced by a transverse field, anisotropy, 
and dimerization in XY spin chains is investigated. The low-energy behavior near the critical point 
is described by a Dirac-type equation for which an exact analytic treatment is possible. Results 
obtained for the dynamical critical exponent, the specific heat, and transverse susceptibihty agree 
with results recently obtained using a real space renormalization group decimation technique. A 
non-zero transverse field changes the universality class of the anisotropy transition. 

Submitted to Physical Review Letters. 
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PACS numbers: 75.10.Jm, 68.35.Rh, 75.10.Nr, 75.50.Ee 

Theoretical studies of quantum phase transitions 
in the presence of quenched (i.e., time- independent) 
disorder has been stimulated by recent experiments 
on spin glasses, ''He absorbed in porous media, and 
superconductor-insulator transitions in dirty thin films. 
Compared to thermal phase transitions in disorder-free 
systems these transitions are poorly understood because 
many of the theoretical methods (e.g., exact solutions, 
renormalisation group, e expansions) that have been so 
useful for pure systems are difficult to implement for dis- 
ordered systems. Yet these phase transitions are asso- 
ciated with particularly rich physics such as large dif- 
ferences between average and typical behavior, new uni- 
versality classes, logarithmic scaling, Griffiths phases j|] 
(in which susceptibilities diverge although there are only 
short-range correlations) , and the breakdown of folklore 
such as "the correlation length is inversely proportional 
to the energy gap" This Letter considers a simple 
exactly-soluble model which has many of these interest- 
ing properties. 

Fisher recently performed an exhaustive study of the 
effect of randomness on the simplest spin model to un- 
dergo a quantum phase transition: the transverse field 
Ising spin chain He used a real space renormaliza- 
tion group decimation technique (RSRGDT) |4| which 
he claims is exact near the critical point. Fisher found 
the phase diagram (which included Griffiths phases near 
the critical point), all the critical exponents, and scaling 
forms for the magnetization and correlation functions in 
an external field. The latter have never been derived for 
the disorder-free case but can be derived in the presence 
of disorder because distributions become extremely broad 
near the critical point. The same model was recently 
studied numerically by Young and Rieger Q , who found 
results consistent with Fisher. The RSRGDT has also 
be used to study the effect of disorder on dimerized |^] 
and anisotropic Q spin chains, chains with random spin 



sizes [||, and quantum Potts and clock chains Possi- 
ble experimental realizations of random spin chains are 
quinolonium(TGNQ)2 and SrgCuPti^^rlrxOg 0- 

This Letter gives an exact treatment of the quantum 
phase transitions in XY spin chains in the presence of dis- 
order. The spin chains are mapped onto a fermion model, 
the continuum limit of which is a Dirac-type equation 
with random mass, for which exact analytic results can 
be derived. The results agree with those of the RSRGDT 
[^,^101, supporting Fisher's claim that it is exact near 
the critical point New results are obtained for the 

anisotropy transition in a non-zero average transverse 
field. It is in a different universality class to the Ising 
transition. 

The Hamiltonian to be considered is: 



i=l 



(1) 



The {erf} are Pauli spin matrices, and the interactions 
Jf , Jf , and transverse fields hi may be independent ran- 
dom variables (see Table I) with a Gaussian distribution. 
The average values will be denoted 



(2) 



and will all be assumed to be positive. 

At zero temperature and in the absence of disorder the 
model undergoes two distinct quantum phase transitions 
jlSf . Both transitions are second order. The transition 
at + jy = h from a paramagnetic to an ferromag- 
netic phase will be referred to as the Ising transition 
The transition at = jy for h < {,r + jv) from 
an Ising ferromagnet with magnetization in the x direc- 
tion to one with magnetization in the y direction will 
be referred to as the anisotropy transition [^^L The 
dimerization transition involves the Hamiltonian (111) with 
(Jf) = (Jf) = J-K (-l)'A and {hi) = 0. At A = there 
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is a transition from a spin liquid to a gapped phase with 
long-range topological order ||^. 

Under a Jordan- Wigner transformation which maps 
spins onto spinless fermions (|l]) becomes 

L 

H= -5]((Jr + Jf)(clc.+i+4iC.) 

i=l 

+ {J^ - Jm\c\^i - c^c^+i) + 1)) (3) 

where boundary terms have been neglected and c| and ci 
denote creation and annihilation operators for a fermion 
on site i. In the disorder-free case (^) can be diagonalised 
by a Bogoliubov transformation ||l3,|l5|. 

The continuum limit of (^) must be taken relative 
to the Fermi wavevector kp^ the wavevector at which 
the energy gap in the fermionic spectrum occurs for 
the disorder-free case This wavevector depends on 
which transition is being considered (see Table I). The 
low-energy properties of all the transitions in Table I are 
described by |16[ 



d 

ivpcr^^ + V(x)a+ + V(x)*a- 
ox 



* (4) 



where cr" (a = 1, 2, 3) and (7± = i(cr^ ±cr^) are Pauli ma- 
trices, ^(x) is a spinor, with (modulo a transformation 
in spinor space) the components given in Table I. V{x) 
is Gaussian white-noise potential with 

{V{x)) = A {V{x)V{x')*) ^-i5{x-x'). (5) 

A measures the deviation from the critical point and |A| 
is the energy gap in the absence of disorder. Table I lists 
values of A, the Fermi velocity vp, and the random vari- 
ables whose standard deviation equals 7, for the Ising, 
anisotropy, and dimerization transitions. For the Ising 
transition (kp = 0) and the anisotropy and dimerization 
transitions in zero transverse field {kp = 7r/2) kp is com- 
mensurate with the lattice and V{x) is real. It will be 
seen that the effect of disorder is significantly different for 
commensurate and incommensurate cases. The contin- 
uum limit describes the low-energy properties [E <^ vp) 
when Vp S> A, ^^7, i.e., arbitrarily close to the critical 
point and for weak disorder. 

It is useful to define an energy D and a dimensionless 
parameter S which are measures of the disorder strength 
and the deviation from criticality, respectively 



Vp 



(6) 



Note that for the Ising transition with — 0, to leading 
order in A/ J^, the parameter S defined by Fisher j|] and 
Young and Rieger ||^ is A/D. 

The advantage of casting the problem in the form of 
the Hamiltonian is that the latter has been stud- 
ied extensively previously, and exact analytic expres- 
sions given for the energy dependence of the disorder- 
averaged density of states {p{E)) and the localization 



length \{E) |T|,|l|,|T|,|^,|l]j2|j2|,||,||. (Due to the 
one-dimensionality all the states are localized by the dis- 
order). The exact results have been found by Fokker- 
Planck equations 18|, supersymmetry 12^,^ , the 
replica trick [ p5[ , S-matrix summation ||20|, and the 
Dyson-Schmidt method . The localization length can 
be found because in one dimension it is related to the 
real part of the one-fermion Green's function [ p2p^ ]. The 
density of states and the localization length can be writ- 
ten in terms of /^(u), the derivative of a dimensionless 
function fs{u) given in Table II, 
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dE \{E) 



inp{E) = npof'siE/D) 



(7) 



where po = l/{nvp) is the value of the density of states 
at high energies {\E\ ^ A,D). 

The low-energy {\E\ <C D) behavior of the density of 
states is given in Table II. For the commensurate case 
the density of states diverges at i? = for (5 < 1/2 and 
is zero at i? = for (5 > 1/2. This difference will lead to 
qualitatively very different behaviour for these two cases. 
In the former case some susceptibilities will diverge as 
the temperature approaches zero. This corresponds to a 
Griffiths phase |l| . The corresponding phase diagram for 
the Ising transition is shown in Fig. ||. The full energy 
dependence of the density of states for various values of 
S is shown in Fig. ||. 

Distinctly different behavior occurs for the incommen- 
surate case. The density of states is always finite at zero 
energy. There is a smooth crossover from gapless behav- 
ior ((p(0)) ^ Po for (5 < 1) to effectively gapped behavior 
((p(0)) <C Po for i5 1). Hence the Griffiths phases still 
exist but no longer have clearly defined boundaries. Sim- 
ilarily, for the commensurate case more than one type of 
disorder (e.g., the anisotropy transition with both a ran- 
dom transverse field and random anisotropy) removes the 
singularity in the density of states . 

The specific heat. Because the eigenstates of the 
Hamiltonian (^ are non-interacting fermions the low- 
temperature behavior of the specific heat and the trans- 
verse susceptibility (for the dimerization and anisotropy 
transitions) follows from the energy dependence of the 
disorder-averaged density of states and are given 

in Table II. The transverse susceptibility clearly shows 
a Griffiths singularity. The results for the commensu- 
rate case agree with results obtained using the RSRGDT 
supporting Fisher's claim |^ that it is exact near 
the critical point. 

Dynamical critical exponent z. This relates the scaling 
of energy (or time) scales to length scales. A crude scaling 
argument |2|] implies that {p{E)) ~ E^/^~^ . Thus for the 
commensurate case, to leading order in (5, z = 1/(2(5), in 
agreement with Fisher |^ and Young and Rieger |^ . This 
is a particularly striking result because it shows that (i) 
z is not universal and (ii) z diverges at the critical point. 
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The latter implies logarithmic scaling and activated dy- 
namics In contrast for the incommensurate case, the 
density of states is finite and constant at low energies and 
so z = 1, as in the absence of disorder. 

Finite size scaling. Monthus et al. [ p9[ studied an equa- 
tion equivalent to (^) with V{x) real and A = [ p5[ . 
They have shown that on a line of length L, for a typ- 
ical potential V{x) the lowest eigenvalue Eq scales like 
-Eq ^ exp(— i^/^). This is consistent with the scaling of 
InE'o with L^/^ at the critical point found numerically by 
Young and Rieger |^. The average (-Eq) ^ exp(— L^/'^) 
p9| , showing the discrepancy between average and typi- 
cal values. 

Correlation lengths. Fisher j|] stressed the dis- 
tinction between average and typical correlations. If 
Cij = (AiAj) denotes a correlation function of a vari- 
able Ai then the average correlation function Cav{r) = 

J2i=i C^i,i+r is what is measured experimentally. Away 
from the critical point Cai,(r) exp(— r/^ai,) where £,av is 
the average correlation length. However, Cav{r) is dom- 
inated by rare pairs of spins with dj ~ 1. In contrast, 
with probability one Ci,i+r ~ exp(— r/^typ) where £,typ 
denotes the typical correlation length. It is distinctly 
different from £^av {£,typ ^ £,av), having a different criti- 
cal exponent. The localization length is useful because 
it is proportional to the typical correlation length for 
quantitites that are diagonal in the fermion represen- 
tation 1^^. Consequently the results in Table II imply 
that for the commensurate case the typical correlation 
length, ^typ ^ '^(O)^^ ^ consistent with previous 

work ]^J31[|. This critical exponent is not modified by 
the presence of disorder. In contrast, for the incommen- 
surate case, A(0) is finite at the critical point. 

Distribution functions. Both Fisher ||^ and Young and 
Rieger ^ considered the distribution functions for vari- 
ous quantities. The present approach can be used to ob- 
tain exact results by using known results for the distribu- 
tion function for the density of states of one-dimensional 
conductors |3|]. The cumulants (or irreducible moments) 
of the transverse susceptibility for the incommensurate 
case at criticality are 



((x..r)c = xo- 



r(|)r(n) 



16T j 



(8) 



where Xo is the susceptibility in the absence of disor- 
der. These moments completely determine the distribu- 
tion function which has been extensively analyzed in Ref. 
^2| . For T ^ D the distribution is Gaussian and cen- 
tered on the mean. In contrast, for T <^ D the distribu- 
tion is extremely broad and asymmetric. The maximum 
of the distribution then occurs near xo exp(— l/s)/2s, 
where s = 16T/(Z?7r), which is much less than the mean 
Xo- This shows the large discrepancy between typical 
and average values. Hopefully, (||) can be generalized to 
the commensurate case and away from criticality. 
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A = h-J''-J'' 

FIG. 1. Phase diagram of the Ising transition in a ran- 
dom transverse field. The horizontal axis is a measure of the 
deviation from criticality in the non-random model. The ver- 
tical axis is the amount of disorder: the rms fluctuation in 
the transverse field. The four phases are: ferromagnet(FM), 
weakly ordered ferromagnet (WO-FM) , weakly ordered para- 
magnet (WO-PM), and paramagnet (PM). The weakly or- 
dered phases are Griffiths phases ||l]. 




FIG. 2. Energy dependence of the disorder-averaged den- 
sity of states for the commensurate case for various values of 
the dimensionless parameter 5 (see Eq. (^)), which is a mea- 
sure of the deviation from criticality. The density of states 
is singular at the Fermi energy [E — 0) when 5 < i. This 
parameter range corresponds to a Griffiths phase. Note that 
only far from criticality {S 3> 1) is there effectively a gap in 
the system. This contrasts with the disorder-free case, for 
which there is always a gap except at the critical point. 
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TABLE I. Different parameters for the continuum limit of three different transitions. Near the critical point all are described 
by the low-energy effective Hamiltonian (|^). The lattice constant is set to unity. The Fermi wavevector kp determines whether 
the problem is commensurate or not. vf is the Fermi velocity. A is a measure of the deviation from criticality. The standard 
deviation of the random variable equals 7 which enters Eq. (p|). 

Transition cos/cf vf A *I'(a; = n)^ Random variable 

lih^i i I J" - .p\ h- -jy (4,c„) h, 

Anisotropy ^ (J" + J^)sinfcF {.F - .P)sinkF (cj, e"'"*^ , c„e-'"'"') Jf - Jl 

Dimerization 2J A {cii-",cii") Ji = Jf 



TABLE IL Summary of the low-energy behavior of different physical quantities. The commensurate case corresponds to the 
Ising transition and to the anisotropy and dimerization transitions in zero transverse field. The incommensurate case corresponds 
to the anisotropy transition in a non-zero transverse field. For the commensurate case qualitatively different behavior occurs at 
the critical point (5 = 0) and away from it {S 7^ 0). The function fs {u) determines the full energy dependence of the density of 
states and the localization length (compare Eq. (W)). Io{x) and Ii{x) are the zeroth and first order modified Bessel functions, 



respectively. r((5) is the gamma function, 
imaginary index. For small x, r{x) ~ x"^ 



/o(a;) 



is a Hankel function of index 5. 
1 + a;^/4, and Ii{x) ~ x. 



Iiy{x) is a modified Bessel function with 
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Incommensurate 
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Localization 
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5— ln(/_2™(25)) 
do 
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